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ABSTRACT 

We present the results from the largest investigation to date of the proximity effect in the H I Lya forest, using 
the fifth Sloan Digital Sky Survey (SDSS) data release. The sample consists of 1733 quasars at redshifts z > 2.3 
showing a signal-to-noise of at least S/N =10. We adopted the flux transmission statistic to infer the evolution 
of the H I effective optical depth in the Lyman forest between 2 < z < 4.5, finding very good agreement with 
measurements from high-resolution quasar samples. By fitting a quasar continuum to each individual object we 
estimated a continuum composite spectrum in the rest frame wavelength range 1000 < A < 3000 A. The flux in 
our composite spectrum closely follows a power law in the form cx v""" with = -0.57 in good agreement 
with similar estimates inferred from significantly smaller SDSS quasar samples. We compared the average 
opacity close to the quasar emission with its expected behavior in the Lya forest and estimated the signature 
of the proximity effect towards individual objects at high significance in about 98% of the quasars. Dividing 
the whole sample of objects in eight subsets according to their emission redshift, we quantify the strength of 
the proximity effect adopting a fiducial value of the cosmic UV background photoionization rate, multiplied 
by a free parameter allowed to vary for different quasars. We inferred the proximity effect strength distribution 
(PESD) on each of the subsets finding in all cases a prominent peak and an extending tail towards values asso- 
ciated to a weak effect. We provide for the first time observational evidence for an evolution in the asymmetry 
of the PESD with redshift. Adopting the modal values of the PESDs as our best and unbiased estimates of 
the UV background photoionization rate (Fhi), we determine the evolution of Fhi within the redshift range 
2.3 < z < 4.6. Our measurements of UV background photoionization rate do not show any significant decline 
towards high redshift as expected by theoretical model predictions and are located at logFei = -11 .78 ± 0.07 
in units of s"' . We employ these estimates and decompose the observed photoionization rate into two major 
contributors: quasar and star-forming galaxies. By modeling the quasar contribution with different luminosity 
functions we estimated their contribution to the cosmic Fhi, allowing us to put a constraint on the remaining 
amount of ionizing photons necessary to obtain the observed photoionization rate. We conclude that indepen- 
dently of the assumed luminosity function, stars are dominating the UV background Fhi at redshift z > 3. 
Subject headings: diffuse radiation - intergalactic medium - quasars: general, absorption lines 



1. INTRODUCTION 

The intergalactic medium (IGM) is permeated by a quasi- 
continuous, spatially fluctuating density field, whose major 
constituents are hydrogen and helium. The distribution of 
these baryons gives rise to a multitude of narrow absorption 
features observed in the spectra of high redshift sources such 
as lum inous quasars (QSOs), coni monly known as the Lya 
forest dRauchll 19981: lMeiksinll2007b . In the framework of hi- 
erarchical structure formation the IGM fills the low density 
regions between collapsed objects and it shows an extremely 
high ionization state. The lack of Lya absorption troughs 
indicates an inhomogeneous distri bution of hydrogen and a 
high ionization state of the IGM dGun n &Peterson|[T965h . 
probed so far out to z ~ 6 dPan et alj 120021) . The cosmic 
UV background (UVB) maintaining this high ionization state 
is generated by the overall population of q uasars and star- 
forming galaxies (Haardt & Madau 1996; F ardal et aljf799l 
iHaardt & Madau .2001 ) . The radiation emitted by flie UV 
sources is filtered by the IGM and diluted by cosmological 
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expansion, leading to a redshift and frequency dependence of 
flie UVB intensity (e.g. lO ave et al.l[T999t iMadau et aT] [T999t 
iSchirber & BtniocHl2003lK 

Three methods are commonly employed to estimate the 
UVB intensity or the UVB hydrogen photoionization rate 
at different epochs. In the first approach, the overall lu- 
minosity density is predicted from a given UV source lu- 
minosity function, taking int o account absorption and re- 
emission by flie IGM (Mada u et al .111 999t IHaardt & Madaul 
120011: ISchirber & Bullock., 2003h . While these computations 
rely on challenging observations such as the galaxy luminos- 
ity function at z > 3, they provide separate estimates on the 
relative contributions of quasars and star-forming galaxies to 
the UVB. 

Secondly, the hydrogen photoionization rate can be in- 
ferred by matching the level of Lya forest abs orption obtained 
in numerical simulations to observations (Rauch et al. 19971 
Theuns et al. 1998; Bolton et al. 2005; Faucher-Giguere et al] 
2008ai) . Although the physics of the IGM is well modeled 
in recent numerical simulations, the primary uncertainty in 
determining the UVB via this method is its degeneracy with 
respect to the observationally poorly constrained thermal state 
of the IGM. 

The arguably most direct observational technique to infer 
the UVB photoionization rate relies on the proximity effect 
of quasars. While fluctuations in the amplitude of the UVB 
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Figure 1. Two examples of quasar spectra taken from the data set with S/N ~ 18 (right-hand panels) and S/N ~ 10 (left-hand panels). In the two upper panels 
the spectra are plotted along with the noise arrays and the power law fit continua. The horizontal arrows over the hya forest identify the wavelength range covered 
in the magnified lower panels where details about the local continua are presented. In all cases the flux density is expressed in units of 10"" erg s"' cm"^ A"' . 
The dashed lines indicate the fitted continua while the soUd Unes show the connected shape as discussed in detail in Sect. |2.2| 



in the general IGM are small (e.g. iMeiksin & Whitjl2004 . 
the UV radiation field near luminous quasars is expected to 
be enhanced. The quasar is capable of over-ionizing the gas 
in its vicinity, resulting in a systematic weakening if the Ly- 
man alpha absorption within several Mpc from the quasar 
dWevmann et al.1 [MTt ICarswell et al.1 [T982t [Murdoch et alj 
I1986I) . The amplitude of this deficit depends on the relative 
contributions of ionizing photons from the local source and 
the UVB. Knowing the QSO luminosit y, the UV background 
photoionization ra te can be estimated (ICarswell et al.lll987t 
lBaithketal.l[T988h . 

Most previous analys es showed a larg e scatt er in the in- 
ferred UVB leve l (e.g. Giallonao et al. 1996; Cooke et aP 
[19971 IScott et al.l [2000; Liske & WilHger 2001), owing to 
the relatively small sample sizes investigated. Although 
the proximity effect is detectable in individu al QSO spectra 
(e.g.[W illiger et al. 1994; Cristiani et al. 1995l lLu et an(l996l; 
ISavadio et al.iil997i;iDairAglio et al.i i2008b). the finite num- 
ber of absorbers in each line of sight gives rise to a substantial 
random scatter that c an only be compen sated with sufficiently 
large samples ( DalFAglio et al.ll2008"al) . 

One of the key questions regarding the UV background 
is its evolution towards high redshift. Observations of the 
proximity effect of a few individual QSOs suggest a mild 
decHne towar ds z ^ 4 (WilH ger et aH 119941; iLu et alJ IT9961 
[Savaglio et aD 119971), in broa d agreement w ith predictions 
from UVB synthesis (Haardt & Madaull2001l) . A similar re- 
sult was obtained by us recently in an investigati on of several 
high-re solution, high S/N quasar spectra (Dall' Aglio et al.l 
l2008al hereafter Paper II). On the other hand, two recent stud- 
ies constraining the UVB from matching simulations to the 
observed Ly forest opacity gave a hydrogen photoionization 
rate nearly constant within 2 < z 5. 4, albeit on rathe r differ- 
ent levels jBolton et al.ll2005l;lFaucher-Giguere et al.ll2008al) . 
It is thus clear that the uncertainties are still substantial. 

In the present study we build on our results obtained in Pa- 



per II where we demonstrated the possibility to use the dis- 
tribution of proximity effect strengths measured in individual 
quasar spectra for constraining the UV background. Our new 
sample of nearly 2000 QSOs from the Sloan Digital Sky Sur- 
vey is by far the largest dataset ever submitted to an analysis of 
the proximity effect, including a significant number of QSOs 
at redshifts z > 4. 

The plan of the paper is as follows. We begin with a descrip- 
tion of the spectroscopic sample of quasars in Section. [H We 
then briefly describe in Section[3]the code adopted to generate 
synthetic spectra. In Section. |4|and|5] we determine the quasar 
composite spectrum and the evolution of the Lya effective op- 
tical depth, respectively. Section |6| introduces the theoretical 
background to measure the proximity effect and presents our 
results on the proximity effect signature towards individual 
objects. In Section [7] we determine the redshift evolution of 
the UVB, followed by a new estimation of the relative share 
of quasars and star-forming galaxies to the UVB. We end with 
our conclusions in Section [9] 

Throughout this Chapter we assume a flat universe with 
Hubble constant Hq = 70 kms"'Mpc~' and density parame- 
ters (f^n^f^A) = (0.3,0.7). 

2. THE QUASAR SAMPLE 
2.1. Selection 

The sample of quasar spectra was drawn from the fifth data 
release ( DRS) of the Sloan Dig ital Sky Survey (SDSS) quasar 
catalog ([Schneider et al.ll2007b . with all data reprocessed by 
the DR6 reduction pipeline. 

The two SDSS double channel spectrographs give a reso- 
lution of FWHM ^2 A(Rr^ 2000) in the wavelength range 
from ~ 3800 A to '--^ 9200 A. We required a coverage of the 
Lyman alpha forest of at least 200 A, corresponding to a min- 
imum emission redshift Zq > 2.3. Thereby we were left with 
~ 18000 of the ~ 77400 visually classified quasars in the 
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Figure 2. Signal-to-noise distribution for each of the quasar subsets as defined in Table[T] 



Table 1 

Subset definitions of the SDSS quasars. 



Subset 


Zq 




Number of quasars 


SI 


2.30-2.50 


2.40 


604 


S2 


2.50 - 2.75 


2.60 


387 


S3 


2.75 - 3.00 


2.88 


153 


S4 


3.00 - 3.25 


3.15 


180 


S5 


3.25-3.50 


3.36 


160 


S6 


3.50-3.75 


3.62 


137 


S7 


3.75-4.00 


3.86 


78 


S8 


4.00-4.611^ 


4.21 


34 



f : The upper limit Zq = 4.61 identifies the highest redshift object 
in this subset 



[Schneider et alJ (12007 1) catalog. In order to be able to reli- 
ably measure the proximity effect in the SDSS spectra, we 
further imposed a minimum signal-to-noise ratio (S/N) of 10 
in the quasar continuum near 1450 A rest frame, reducing the 
number of QSOs to 1916. Among these there are 144 broad 
absorption line quasars (BAL QSOs) which cannot be used 
for the proximity effect ana lysis. Details about their identifi- 
cation are given in Sect. 12. 31 Finally we inspected the remain- 
ing quasar spectra visually and removed 39 objects that were 
misclassified or were still suspected to be BAL QSOs. This 
procedure left us with a total of 1733 quasars. 

The proximity effect analysis also requires accurate quasar 
redshifts and s pectro photometry. The redshifts in the 
[Schneider et alJ (l2007h catalog were measured from broad 
emission lines in the accessible wavelength range and sta- 
tistically corrected for the systematic velocity shifts between 
low- and high-ionization lines (e.g. Gaskell 1982). The sta- 
tistical correction was determined bv iShen et al.l (i2007 ) from 
a quasar sample at 1.8 < Zq < 2.2, where the SDSS spectro- 
graph covers the Si IV,C IV, C III], and Mg II emission lines. 
Adopting the Mg II emission redshift as systemic and apply- 
ing their velocity correction, the redshift accuracy of higher 
redshift quasars is still azq — 0.002, sufficient for our pur- 
poses. Spectrophotometric accuracy is achieved by observing 
calibration stars together with the SDSS survey targets and ty- 
ing the spectra to the accurate SDSS photometry (crmag = 0.04; 



lYork et alJl2000l) . 

In the following we divide the sample into 8 subsets accord- 
ing to the quasar emission redshift as listed in Table [T] Each 
subset is characterized not only by the range in redshifts, but 
also its number-weighted mean Zmean- As an example of the 
average data quality, we show in Figure [T] two quasar spec- 
tra together with their la noise arrays and the estimated con- 
tinua which will be discussed in the next section. Addition- 
ally, Fig.|2]shows the S/N distribution for the subsets. 

2.2. Continuum estimation 

In order to evaluate the proximity effect along each sight 
line, the quasar flux at the Lyman limit and the shape of its 
continuum emission are needed. Therefore two types of con- 
tinua were estimated: (i) a global power law {fu{v) oc t^"), 
excluding emission and absorption regions, used to esti- 
mate the quasar flux at the Lyman limit and (ii) a more lo- 
cal estimate that also includes the broad emission lines as 
a quasi-continuum, used to infer the continuum-normalized 
absorption-line spectrum. 

The power law continuum was fitted though pre-defined 
continuum regions redward of the Lya emission line free of 
any known emission lines (Sect. [Hi. Extrapolating the power 
law to the Lyman limit, we estimated the quasar Lyman limit 
flux. The necessity for this extrapolation is motivated by the 
difficulty to infer the Lyman limit flux due to the low S/N at 
short wavelengths or the lack of wavelength coverage at low 
redshifts. 

The local continuum was auto matically fitt e d usin g an 
algorithm based on the work by lYoung et alj ( 11979 ) and 
ICars well et alJ ( Il982h . A cubic spline was interpolated on 
adaptive intervals along the spectrum with respect to the con- 
tinuum slope. The points for the spline interpolation were 
chosen starting from a regular sampling of the spectrum with 
a binning that becomes finer whenever the slope of the com- 
puted continuum exceeds a given threshold. This was done in 
order to better reproduce the wings of emission lines. To this 
continuum we additionally applied a systematic correction 
which accounts for the resolution effects and line blending, 
estimated from Monte-Carlo simulated spectra (see Sect. [3] 
for details). We basically followed th e procedure outlined in 
Paper II and lDaU'Adio etalj ( l2008bl) hereafter Paper I. This 
approach consists of analyzing the ratio between the input and 
fitted continuum in a set of 500 mock spectra. From this ratio 
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Figure 3. Average ratios between the fitted and input continua in the Lya 
forest of our simulations (red lines). The smooth profile with the dashed lines 
represents the least-square second order polynomial fit to the ratios and the 
± 1 (T eiTors. 

we determined an average continuum ratio and a dispersion 
used as systematic and statistical uncertainties. Our results 
are presented in Fig. [3] The local, corrected continuum was 
finally used to estimate the quasar transmission T = /\qso/^cont> 
defined as the ratio of the quasar flux (Fqso) and the local con- 
tinuum (Fcont)- 

2.3. Broad absorption line quasars 

The input sample of QSOs contains a significant fraction 
of BAL quasars which are not suitable to study the proximity 
effect and therefore need to be removed. To i dentify these ob- 
j ects w e followed the procedure outlined by IWevmann et alj 
(119911) . The approach consists of computing the balnicity in- 
dex defined as 

BI = - I [l-r(v)/0.9]Cdv (1) 

with vi = 25000 km s"^ and V2 = 3000 km s"' . T{v) is the trans- 
mission as function of velocity and C is a constant initially set 
to zero. C becomes 1 .0 whenever the term in the brackets re- 
mains positive for at least 2000 km s"' , meaning that the first 
2000 km s"' of a BAL do not contribute to the balnicity index. 
C returns to zero whenever the quantity in brackets becomes 
negative. If a quasar shows a positive balnicity index, it was 
flagged and removed from our sample. 

We note that in spite of the seemingly objective definition 
of the balnicity index, the identification of a BAL quasar still 
depends on the subjective choice of a number of factors such 
as the quasar continuum or the velocities vi and V2. With this 
procedure we identified 144 BAL quasars not suitable for the 
analysis of the proximity effect. We noticed later that a few 
objects with BAL-like features survived this deselection and 
remained in the sample . We discuss their role in the proximity 
effect detection in Sect l6.3l 

3. MONTE CARLO SIMULATIONS 

Simulating quasar spectra provides a powerful tool of com- 
parison which we employ to accurately control and calibrate a 
number of systematic effects present in real observations. We 
employed such simulations for three main purposes, (i) We 
estimated the systematic and statistical uncertaint ies o f our 
continuum placement in the quasar spectra (Sect. I2.2l i. (ii) 
These uncertainties were used to obtain error bars on the ef- 
fective optical depth in the forest, in addition to the variations 
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Figure 4. The flux-median and continuum-mean SDSS composite spectra 
derived from our set of QSOs. The vertical dotted lines mark the positions of 
the most prominent emission lines known in quasar spectra. The difference 
between the two composites blueward of the Lya emission line is due to the 
influence of the Lya forest absorption. 

in Lya absorption among different lines of sight (Sect. |5]). 
(iii) We used the simulated spectra to assess the bias aris- 
ing from an asymmetric proximity effect strength distribution 
when combining the data to estimate the UVB photoioniza- 
tion rate (Sect.|71l. 

The procedure used to generate synthetic spectra is identi- 
cal to the one presented in Paper II. We generated a set of 500 
synthetic spectra adopting a simple Monte Carlo prescription, 
assuming that the Lya forest is well represented by three ob- 
served distributions: 

1. The line number density distribution, approximated 
by a power law of the form dn/dz(x{\+zy where 
7 = 2.04 (Paper II). 

2. The column density distribution, given by 
/(Nm ) c^N^^^ where the slope is /3 ~ 1.5 CKim et al.1 
l200ll) . 

3. The Doppler parameter distribution, given by 
dn/db oc b~^sxp \-bt/b'*] where ^j^r— 24km/s 
(iKim et al.ll200lh . 

The algorithm populated the simulated sight lines with ab- 
sorption features following the above distributions until the 
evolution of the effective optical depth matches the estimates 
made in Paper II. Once the simulated transmission was com- 
puted, an artificial quasar spectral energy distribution includ- 
ing emission lines of varying strengths and widths was gener- 
ated via the principal component method described by Suzukil 
(2006) and in Paper I. The resolution of the spectrum was de- 
graded using a Gaussian smoothing function with a FWHM of 
approximately 2A and pixel size increased to match the SDSS 
specifications. Gaussian noise was added to the final quasar 
spectrum, in order to reproduce the S/N level of our objects. 

4. A QUASAR COMPOSITE SPECTRUM 

The large number of objects in our sample motivated us to 
compute a composite quasar spectrum. This spectrum pro- 
vides important insights about the position of low and high- 
ionization lines as well as global properties of the quasar spec- 
tral energy distribution. The computation of the composite 
spectrum required: (i) A rebinning of the spectra in rest frame 
wavelengths, (ii) a normalization of the fluxes and (iii) a com- 
bining method to obtain the final composite spectrum. We 
computed two types of composite spectra, each of which has a 
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Figure 5. The dashed Hne represents the power law fit {fu(y) oc i^™) to the 
continuum-mean SDSS composite spectrum. 



different combining procedure: (1) the median spectrum em- 
ploying the fluxes of the QSOs to minimize the effect of ab- 
sorption Unes and (2) the mean spectrum adopting the fitted 
continua of the QSOs. In the latter approach all quasar con- 
tinua were corrected for th e sys tematic depression in the Lya 
forest as described in Sect. 12.21 Both results are presented in 
Fig- SI The two composites are essentially identical until the 
Lya emission where the differences between the two become 
visible. These differences are due to the absorption in the Ly 
forest which we accounted for when combining the corrected 
quasar continua. 

The very high S/N level of our continuum composite spec- 
trum allowed us to identify many of the well known emis- 
sion lines present in the spectra of high redshift QSOs. To 
confirm the accuracy o f the quasar redshifts provided in the 
[Schneider et alj (120071) catalog, we measured from the com- 
posite spectrum the location of the most prominent low and 
high-ionization emission lines. For lines such as the Lya, 
Si IIH-O I, Si IVH-O IV], C IV, C III, and Mg II we fitted one or 
more Gaussian profiles to accurately estimate the emission 
wavelengths. The measured values for low-ionization lines 
were in agreement with their laboratory wavelengths while 
the high-ionization lines present ed a systematic blueshift up 
to 550 km/s in flie case of C iv (lGaskelllll982l: iTvtler & Faij 
[1991 [Richards et al.ll2006h . 

The computation of the two composite spectra enabled us to 
check the accuracy of the BAL quasars removal. If a signif- 
icant fraction of such objects were misidentified, we should 
have detected a depression only in the flux composite, red- 
ward of the C IV emission line. The lack of such a feature 
ensures a correct removal of any significant contamination by 
BAL quasars. 

The mean continuum composite spectrum is plotted in 
Fig. H] with a single power law fit. The accuracy of this 
fit is compromised by the limited number of regions in a 
QSO spectrum free of emission lines. To perform the power 
law fit we adopted the three regions between 1350-1365 A 
1445-1475 A and 2200-2250 A. The estimated spectral slope 
is = -0.57, in agreemen t with previous measure ments 
based on earlier data releases jVanden Berk et all20011) . This 
value is further supported by the distribution of spectral in- 
dices which closely follows a Gaussian distribution centered 
at = -0.6 with a dispersion of (Ja = 0.35. Therefore, esti- 
mating the quasar flux at the Lyman limit via a single power 
law extrapolation ensures reasonable accuracy in the compu- 
tation of the proximity effect. 



Table 2 

Effective optical depth in the Lya forest of our combined QSO sample, 
connected for the continuum uncertainties. 



< z > 


f"cff 




< z > 


Teff 




2.16 


0.183 


0.002 


3.36 


0.531 


0.009 


2.26 


0.211 


0.003 


3.46 


0.561 


0.010 


2.36 


0.233 


0.003 


3.56 


0.617 


0.011 


2.46 


0.256 


0.003 


3.66 


0.668 


0.012 


2.56 


0.291 


0.004 


3.76 


0.729 


0.016 


2.66 


0.315 


0.004 


3.86 


0.748 


0.016 


2.76 


0.331 


0.004 


3.96 


0.804 


0.017 


2.86 


0.359 


0.008 


4.06 


0.787 


0.025 


2.96 


0.376 


0.008 


4.16 


0.888 


0.026 


3.06 


0.426 


0.008 


4.26 


0.895 


0.027 


3.16 


0.459 


0.008 


4.36 


1.043 


0.028 


3.26 


0.499 


0.008 









5. EFFECTIVE OPTICAL DEPTH OF THE LYa FOREST 

The signature of the proximity effect is typically observed 
against the increase with redshift of the mean absorption in 
the Lya forest. The opacity of the Lya forest is empirically 
known to scale with redshift according to a power law of the 
form 

Teff=ro(l+zr^'. (2) 

for a line density evolving as dn/dz oc (1 -\-zV ■ In Equation|2] 
Teff is the H I effective optical depth. This quantity is related 
to the mean transmission as Teff = -In (T), where the average 
( ) is taken in redshift intervals. 

The adopted method to estimate the parameters (logT(),7) 
consists of measuring Teff in the forest of those quasars which 
intersect a particular redshift slice (Az = 0.1). We removed 
the influence of the quasar by considering only the rest 
frame wavelength range 1025 < A < 1180 A. We performed 
this computation first on the transmission without continuum 
correction, and then applied the correction as described in 
Sect. 122] 

The results are shown in Fig. |6] and listed in Table |2l It 
is evident that the continuum correction is crucial to prevent 
strong biases in measuringthe evolution of Teff at low reso- 
lution. The best fit of Eq.^to the continuum-corrected data 
yields values of (logTo,7) = (-2.28 ± 0.03,2.13 ± 0.04). De- 
spite the moderate spectral resolution, our results are in ex- 
cellent agreement with the evolution of the Lya optical depth 
as recently measured by us from high resolution quasar spec- 
tra in Paper II and also by different authors (e.g. lSchave et alj 
120031: lKimetani2007l) . We also observe that for z > 3.5 our 
estimated Teff(z) is significantly lower t han recent measure- 
ments by 'Faucher-Giguere et alj (l2008bl) . based on a sample 
of 84 QSO spectra obtained by different instruments. We have 
no explanation for this discrepancy. Because of the excellent 
agreement with other studies we proceed adopting our newly 
estimated fit parameters for further computations. 

As a parenthetical note we remark that we do not detect 
any deviation from a single power law evo lution at z ~ 3.2. 
Around this redshift, iBernardi et alj (I2003h detected a dip in 
7"eff(z) based on a much smaller sample of SDSS quasar spec- 
tra and speculated that this might be related to He II reioniza- 
tion. While our higher resolution data do show a similar dip 
(cf. Paper II), there is no evidence for such a feature in Fig.|6] 
A more detailed analysis of the reality of that dip is probably 
worthwhile, but outside the scope of the present paper. 
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Figure 6. Observed effective optical deptli as a function of redsliift witli and without the continuum correction (dots and triangles, respectively). The solid 
profile shows the power-law least-squares fits to the data points while the dashed, dotted and dashed-dotted curve indicates the evolution of the hya optical depth 
measured in Paper II and, by Faucher-Giguere et al. (2008b) and Kim et al. 12007) respectively, using high resolution spectra. Errors on the redshift axis indicate 
the portion of forest used to compute the effective optical depth. The uncertainties on the effective optical depth have been estimated using our Monte Carlo 
simulated sight lines. 



6. MEASUREMENTS OF THE PROXIMITY EFFECT 
6.1. Definitions and procedure 

The ionization state of the IGM near a quasar will be dom- 
inated by the local source rather than the UVB, leading to 
the observed weakening of the Lya absorption in its vicinity. 
An efficient technique of detecting the proximity effect is to 
compare the effective optical depth measured nea r a quasar to 
its expected value in the Lyman alpha forest (e.g. lLiske et alj 
[1998; Liske & WilHger 2001). The flux statistic technique is 
applicable to high- and low-resolution spectra and has been 
used by us to detect the proximity effect on individual lines of 
sight (Paper I & Paper II). 

The influence of the quasar modifies the effective optical 
depth to 

T,ff = To(l+zr^\l+co)'-f^ (3) 

where /3 is the slope in the column density distribution. The 
quantity ui is defined as the ratio between the quasar and back- 
ground photoionization rates: 



^(z) = 



1 



rHi(z) rH:(z) 



hv 



a{v)dv. (4) 



We assume the UVB photoionization rate to be reasonably 
constant for the redshift path length probed by the Lya forest 
of a given QSO. Approximating the spectral energy distribu- 
tion of a quasar in the form fuW) oc we can rewrite Fq 
as: 

y , gz.(i^Q) L^jvo) ( 1+Zq \ """" 

"f^^^ A^hO-a^)dl(z^,z)\l+z) ^' 

where z is the redshift along the line of sight such that z < Zq, 
dL(Zq,z) is the luminosity distance of the QSO as seen from 
any foreground redshift along the LOS, is the spectral 
slope of the QSO and L^^, its luminosity at the Lyman limit. 
The hydrogen photoionization cross section at the Lyman 
limit is (Tuivo) = 6.33 x 10"'^ cm^ and h is the Planck's con- 
stant. 

We emphasize that all previous analyses of the proximity 
effect were based on a simplified determination of u)(z) as a 
flux ratio by assuming the same spectral energy distribution 
of the quasar and the cosmic UV background in Equation |4] 
With our correct definition of omega as a ratio of photoioniza- 
tion rates, we do not have to assume a UVB spectral energy 



distribution to get poorly constrained UVB intensities. More- 
over, all recent numerical studies of the UVB give the UVB 
photoionization rate, facilitating comparisons between differ- 
ent methods. 

Finally, the ratio of the observed optical depth to the one 
expected in the Lya forest, or the normalized effective optical 
depth ^, is given by 



1-/3 



(6) 



ro(l+zr+i 

with the parameters (logro,7) measured in Sect. |5] The 
slope of the column density distribution was fixed to /3 = 1 .5 
(Kim et al. 200li). 

The proximity effect will manifest its signature as a depar- 
ture from unity in the normalized optical depth as cj ^ cxd, i.e. 
towards the quasar systemic redshift. 

6.2. The proximity effect along single sight lines 

We now explore the application of the above method to in- 
dividual objects. Our approach is as follows: We first com- 
pute, for each individual line of sight, the normalized effec- 
tive optical depth as a function of lo, within bins of logcj. We 
then search for a systematic decrease of ^ at high values of uj. 

In order to compute the uj scale (Eq. Hjl, a fiducial value for 
the UVB photoionization rate had to be assumed. We fixed 
this value to F* = 10"'^ s"', thereby uniquely converting the 
redshift scale into an u) scale. Dividing the logcj axis into an 
equally spaced grid, we then modeled the decrease of ^(uj) in 
bins of Alogw according to the formula 



1-/3 



(7) 



where a is the only free parameter which expresses the ob- 
served turnover of ^ in units of F*. In principle, each a di- 
rectly provides a best-fit value of Fhi, since Fhi = a x F*. As 
demonstrated in Paper II, however, the random scatter of a 
due to the finite small number of absorbers in each line of 
sight is so large that meaningful constraints on Fhi can only 
be obtained from the distribution of a values for a statistically 
significant sample. 

For a given sight line, the value of the fitting parameter a 
describes how much the best-fit model disagrees with the fidu- 
cial value assumed above. We regard logo as a measure of the 
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Figure 7. Thi'ee examples of the proximity effect signature along individual sight lines at redshift z ~ 3. Each panel shows the normalized effective optical depth 
5 versus uj, binned in steps of A log a; = 1, with the fiducial model superimposed as the dotted line. The solid lines delineate the best fit to each individual QSO 
as described in the text. This subset was chosen for presentation purposes to show the variable strength of the proximity effect, going from strong (left panel) to 
weak (right panel). 



Strength of the proximity effect, which can be classified into 
three principal regimes as presented in Fig.|2l 

1 . log fl < 0: this case describes a strong proximity effect. 
The normalized effective optical depth is already small 
forw < 1. 

2. logfl ~ 0: this case describes an average proximity ef- 
fect since the ^ values follow the fiducial model profile. 

3. log a > 0: this case describes a weak proximity effect. 
The normalized effective optical depth will turn over 
only at very large logw values (if at all). 

Of course, the boundaries between these three categories de- 
pend on the (arbitrary) choice of the fiducial UVB photoion- 
ization rate value F* which fixes the zero point of log a. 

6.3. The distribution of proximity effect strengths 

We quantified the proximity effect in the above manner on 
each sight line in the sample to compute the distribution of 
logo. Considering the UVB as a constant for each subset of 
Table [T] we plot the observed proximity effect strength distri- 
bution (PESD) of each subset in Fig.H] 

Each of these distributions shows a wide spread in log a val- 
ues. For in total 35 objects, or 1.8 % of the sample, we obtain 
logfl > 3.7, which means that for these objects the fit hinges 
on a single point showing a decrease of £^{uj) at large uj. Since 
the largest values of log a correspond to the weakest proxim- 
ity effect signatures, we inspected these cases individually. It 
turned out that all of them displayed either strong associated 
absorption or were previously unrecognized BAL QSOs. For 
all other QSOs in the sample, the proximity effect signature is 
based on at least two and in general several points in the ^ vs. 
logo; diagram; we conclude that the effect is detected in 98 % 
of the sample. 

For all subsamples the PESD shows the two major char- 
acteristics discovered in Paper II: (i) A well defined peak, 
and (ii) a significant asymmetry with an extended tail towards 
weaker proximity effects, i.e. large values of a. Additionally, 
the PESD becomes substantially more narrow with increasing 
redshift and also shows a less prominent tail towards log a > 0. 

The skewness in the PESD is the result of the uj definition as 
a function of redshift (Eq.H]). Approaching the quasar emis- 
sion, equal A log uj bins probe progressively smaller Az inter- 
vals, leading to a deviation in the distribution of the absorber 
counts from a Gaussian. This translates to a skewed distri- 
bution of ^ as z Zq since large ^ values (and therefore a 
weaker observed proximity effect) become more likely. This 



trend becomes less prominent at high redshift because equal 
Alogix) bins will probe progressively increasing numbers of 
absorbers and therefore reduce the skewness in the PESD. 

The PESD asymmetry is further amplified by the presence 
of overdensities near some of the QSOs in each subset. In Pa- 
per II we found that about 10% of our objects showed a signif- 
icant excess of absorption within a few Mpc of the quasar red- 
shift. The inferred values of a for these quasars are expected 
to be larger by a factor related to the amount of absorption 
which becomes observable in the PESD as an enhancement 
in the tail extension on its weak side. Since we do not have 
sufficiently high resolution in our spectra to estimate the ex- 
cess of absorption on Mpc scales near the quasar emission, we 
estimated its influence by employing our Monte Carlo simu- 
lation. 

We generated a set of synthetic quasar spectra as described 
in Sect. [3] matching the characteristics of each of our sub- 
sets (sample size, and redshift and luminosity distributions). 
For each spectrum we performed a continuum fit, also ap- 
plying the appropriate correction and then we estimated the 
proximity effect strength. This exercise serves two purposes: 

(i) as an accuracy check for our automatic data analysis, and 

(ii) as a tool to reveal the effect of H I overdensities on the 
PESD. For the input photoionization rates in the different sub- 
sets we a dopted the modal value of the observed PESDs (see 
Sect.|73]l. 

Each panel of Figure [8] shows the comparison between the 
observed and the simulated PESD. These distributions clearly 
differ with respect to the tail extension towards weaker prox- 
imity effect, while they are in agreement with respect to the 
location of the most likely value of logo. The simulated PESD 
is narrower than the observed one, indicating that a certain 
fraction of QSOs may be affected by large-scale overdensities. 
Our approach to estimate the UV background from these dis- 
tributions rests on the assumption that the peak of the PESD 
is unaffected by this contamination. We discuss the validity 
of this assumption further below. 

7. THE EVOLUTION OF THE UV BACKGROUND 
PHOTOIONIZATION RATE 

7.1. Method 

As shown in Paper II, the shape of the PESD depends on the 
cosmic photoionization rate as well as on possible excess H I 
absorption on scales of several Mpc around the quasars. We 
now estimate the photoionization rate Fhi in various redshift 
bins by evaluating the PESD. 

The first method consists of recovering the modal value of 
the PESD, i.e. estimating the peak of the distribution. To 
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Figure 8. The proximity effect strength! distributions in individual subsets of our data set. Each panel shows the observed PESD for increasing redshift intervals 
(black histogram) along with the simulated PESD obtained from a set of simulated lines of sight (gray histogram). The dotted line and the point mark the modal 
values of the PESD with the associated uncertainties. 



achieve this goal we adopted the following bootstrap tech- 
nique: Starting from the observed distribution of A^, values of 
logo, where A^, is the total number of objects in a given sub- 
set, we randomly duplicated the strength parameter of Ni/e 
quasars and estimated the modal value of the new PESD. We 
repeated this process 500 times in each of the eight subsets, 
obtaining the mean and the sigma values of PESD modes. 

In addition to the modal values of the PESD, we also ex- 
plored another method to estimate the photoionization rate 
evolution with redshift. We take advantage of the fact that 
overdense quasars will have statistically larger logo values 
and selected from our sample only those QSOs with logo 
below a certain threshold logoo- We fixed logao =1.5 af- 
ter comparing how the simulated and observed PESD change 
with redshift (see Fig.[8]l. This is equivalent to truncating the 
PESD above logao- Such a cut does not modify the estimates 
of the PESD mode, but significantly reduces the biases con- 
nected to overdense quasars. We then combined the signal of 
all these quasars by merging them according to their redshift 
scale. We derived the ^ values in the same way as for mea- 
suring the proximity effect along single sight lines and then 
fitted Eq. |2]to the data, yielding a first estimate of Fhi. We 
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Figure 9. Comparison of the two different estimates of the UVB photoion- 
ization rate inferred from the modal value of the PESDs (solid dots) and 
from the combined analysis of the proximity effect after truncating the PESD 
as described in the text (empty squares). We also show the best-fit non- 
evolutionary model to our data (dashed line) yielding a photoionization rate 
of log Thi = -1 1 .78 ± 0.07 in units of s"' . 



finally corrected the preliminary photoionization rates for the 
systematic biases introduced by the remaining PESD asym- 
metry. 

Figure|9]shows a comparison between our two estimates of 
Fhi- The UVB values derived via both methods are consistent 
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Figure 10. Evolution of the UVB photoionization rate as obtained employing the modal value of the PESDs on the SDSS data set (solid circles ) . The 
empty squares refer to the estimates from the UVES sample presented in Paper II. We compare our results with recent estimates from I Bolton et al.l 120051) 
and Faucher-Gieuere et al. 1 2008a), both employing numerical simulations to infer the cosmic photoionization rate (crosses and empty circles, respectively). 



with another, yielding a roughly constant level of the H I pho- 
toionization rate at 2 < z < 4.6. A least-square fit to a non- 
evolving photoionization rate yields log Fhi = -11 .78 ± 0.07 
for the PESO modal values and log Thi = -11 .75 ± 0.05 from 
the second method, both in units of s"'. While the two meth- 
ods yield statistically very similar results, we adopted the 
modal estimates of the UVB photoionization rate as our fi- 
nal estimates. This choice is mainly motivated by the concep- 
tual simplicity of the PESD mode and its independence of ad 
hoc assumptions regarding the threshold logoo or the PESD 
asymmetry correction. 

7.2. Effect of overdensities on the proximity effect strength 
distribution 

We now address the concern that large-scale H I overden- 
sities around our quasars might significantly bias the UVB 
estimates. A fixed overdensity factor valid for all quasars in a 
sample would simply shift the PESD towards larger log a val- 
ues, but preserve the shape of the distribution. Our method is 
completely degenerate to such a shift. We argued in Paper II 
that a uniform overdensity factor is however rather unplau- 
sible. In the more likely case of a distribution of overdensi- 
ties the width of the PESD would be significantly enhanced. 
Measuring the width of the distribution thus provides a useful 
additional diagnostic. 

In Fig.[8]we compare the measured PESDs with simulated 
distributions generated by the Monte-Carlo process as de- 
scribed above. In each panel, the simulated dataset has the 
same redshift distribution as the observed one, and the areas 
under the observed and simulated histograms are normalized 
to the same total number of objects. 

In most of the panels of Fig. |8] the two histograms agree 
quite well in the cores, but the observed data show a signifi- 
cant excess of large logo values that is absent in the simula- 
tions. In terms of standard deviations, the observed distribu- 
tions have dispersions typically ^ 1 .5 x that of the simulated 
ones; most of that excess is due to the objects with log a > 2. 
We interpret this excess as the main manifestation of H I over- 



densities around our quasars. 

The applicability of our approach to estimate the UVB via 
the proximity effect now rests on the question whether or not 
the mode of the PESD is significantly affected by the influ- 
ence of overdensities. The generally good agreement between 
the cores of the simulated and observed distributions suggests 
that only a minority of QSOs produce the excess; this supports 
our notion that the mode is approximately conserved. 

Consider our highest redshift subset S8 in particular, con- 
taining 34 QSOs at z > 4. The PESD of this subset has a 
dispersion of a\oga = 0.51, while the corresponding simulated 
PESD has (T = 0.46. Essentially all of this difference is due to 
a single outlier. Any excess H I absorption in the QSOs would 
make the PESD broader and inconsistent with the observed 
narrow distribution, unless the amount of such an excess ab- 
sorption were nearly exactly equal for the entire sample. 

We conclude that the PESD mode is quite robust against 
the influence of overdensities. While we cannot firmly ex- 
clude the possibility that some of our Fhi values at interme- 
diate redshifts are somewhat affected, the narrowness of the 
z > 4 PESD provides evidence against strong H I overdensi- 
ties on scales of several Mpc. We consider this data point as 
our most significant measurement. 

7.3. Comparison with other measurements 

FigurefTOlshows our estimated values of the hydrogen pho- 
toionization rate Fhi as a function of redshift, together with 
recent results from the literature. In order to facilitate a di- 
rect comparison, we repeated the proximity effect analysis of 
our high-resolution VLT/UVES spectra (Paper II), but now 
quoting photoionization rates rather than UVB intensities and 
applying a coarser binning. The redshift-averaged Fhi values 
are very similar, -11.78 ±0.07 for SDSS vs. -11.83 ±0.10 
for UVES (in logarithmic units of s"'). While at z < 3 the 
two datasets are almost indistinguishable, the UVES measure- 
ments suggest a decrease of Fhi with increasing z while the 
SDSS data appear to favor an approximately constant pho- 
toionization rate. The differences are not highly significant 
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considering the error bars; we also recall that the UVES sam- 
ple contains only 7 QSOs with z > 3.3. 

The SDSS measurements are also in broad agreement with 
the results by Bolton et al. ( 2005) who inferred Fhi by match- 
ing hydrodynamical simulations to the observed opacity evo- 
lution in 19 QSOs (with only 4 located at z > 3.3). Bolton 
et al. discussed the main uncertainties entering in the compu- 
tation such as the IGM temperature, the IGM density distri- 
bution and the evolution of the effective optical depth. Their 
best values of Fhi are slightly smaller than our SDSS esti- 
mates, but again consistent within the error bars. In particular, 
their results also show no downturn of the photoionization rate 
towards large redshifts. On the other hand there is less agree- 
me nt between our measurements and the latest estimates by 
by iFaucher-Giguere et al.l ( l2008ah . which are globally offset 
by ~ 0.5 dex towards low values with respect to the SDSS 
results (but again show no trace of a downturn). This discrep- 
ancy is presumably (at least partly) due to the already men- 
tioned discrepant description of the optical depth evolution 
Teff(z) (see Sect.llll. 

These estimates of the UVB can be compared to the semi- 
analytic prediction based on UV luminosity functions of 
quasars an d star-for ming gala xies ( Haardt & Madau 1996; 
iFardal et al . 1998; H aardt & M adau 200JD. As presented in 
Fig.[lOl the theoretical calculations done bv lHaardt & Madaul 
(HOOl) predict a significant decline of Fhi towards high red- 
shifts, being consistent within the errors with both the results 
of Bolton et al. (2005) and our own from UVES. Since empir- 
ical constraints on both the quasar luminosity function and on 
the high-redshift stellar emissivity have been much improved 
in recent years, we present in the following a new calculation 
of the UVB. 

8. ORIGIN OF THE METAGALACTIC UV PHOTONS 
8.1. Method 

The intensity of the background radiation field at a given 
frequency and redshift is derived from the cosmological ra- 
diative transfer equation 



d_ 
di 



vH- 



d_ 

dv 



Ji, = -3HJ^ -cavJ^ + —e^ 

4-TT 



(8) 



( lPeebleslll993h . where H(t) is the Hubble parameter, c is the 
speed of light, and and represent the absorption and 
emission coefficients, respectively (note that a^, has no con- 
nection with the quasar slope estimated in Sect.Hji. The aver- 
age intensity of the radiation can be expressed, at an observed 
frequency i^o and at a given redshift zo, as 



A-TT 



e.(z.,z)e-^-(''°'^°'~->^dz (9) 



dz 



with I' = i'o(l +z)/(l +Zo), dl/dz being the proper length ele- 
ment and TLyc is the Lyman continuum opacity, describing the 
radiation filtering due to intervening absorption. 

The rapid increase of opacity along the LOS for z > 2 
leads to a significant reduction in the mean free path of ioniz- 
ing photons, which as a consequence implies that only local 
sources contribute significantly to the UV background inten- 
sity (Madau et al. 1999; Schirber & Bullock 2003). In this ap- 
proximation, the solution of Eq.|9]is given by 
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Figure 11. QSO emissivity ev(z) at the Lyman limit obtained from four 
different luminosity functions as described in the text. The different lines 
refers to [a] the Bongiorno et al. 1 2007) LDDE model, [b] the modified LDDE 
model, [c] and [d] the iWolf et aU 12003.) PLE and PDE model, respectively. 

The mean free path /S.I{vq,z) is therefore related to the prop- 
erties of the absorbers in the Lya forest at redshift z > 2 and 
is given at the Lyman limit frequency by 



1 



(11) 



where F(2 - /3) in this case is the Gamma function and 
(Nq,/3,j) describe the redshift and column density distribu- 
tions of the absorbers, typically well represented by 



dz dNn, 



(12) 



The major contribution to the mean opacity at the Lyman 
limit is driven by systems with optical depths around unity, 
or equivalently with column densities of the order of A^hi — 
10'^'^ cm"^. In this regime the properties of the absorbers 
in the Lya forest are particularly diffi cult to infer and only 
few l i terature measurernents exist (e.g.lStengler-Larreaet alj 
IT99I iKiAinan & Tvtlerill997l; iPeroux et al.ll2003l) . For sim- 
plicity and due to the non-negligible contribution of lower 
column density systems to the Lyman continuum opacity 
(iMa dau et aljl99 9). we approximated the column density dis- 
tribution to a single power law with {3 = 1.5. The absorber 
density evolution has been fixed to a power law as well with 
(Na.-y) = (1.4 X 10^2.45). following |P6roux et alj (120031). De- 
viations from these assumptions will only have a marginal 
influence on both the global and relative scaling of different 
contributions which, as we will show, does not affect our re- 
sults. 

Assuming that the UV background intensity has a power- 
law shape for wavelengths shorter than the Lyman limit (Ji, oc 
I/"''), the UVB photoionization rate becomes 



rHi(z): 



a-Aj^o) A/(i/o,z) e^.(^'o,z) 
h(3- ah) 



(13) 



The spectral index of the UV background, a^, is a poorly 
known quantity which in principle could be predicted as- 
suming a s pecified contribution from quasars and galaxies 
(iHaardt & M adau 20 03. In our calculation we fixed its value 
to ab = -1 .4 ( Ag afonova et al.ir2005l) . but note that a different 
value would not change our conclusions as will be discussed 
in Sect. 18.41 Equation[T3]directly converts the integrated lumi- 
nosity density at the Lyman limit, e^i^jZ), into a photoioniza- 
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Figure 12. QSO contributions to the UVB photoionization rate (labeled 
lines) as obtained from the dilferent emissivities presented in Fig. llll in com- 
parison to our estimates of the total photoionization rate from the SDSS (solid 
circles) and the UVES samples (empty squares). The solid line represent the 
assumed level for the total photoionization rate of IorFhi = —11.7 8. The 
different curves were inferred form [a] the Bongiomo et alj 120071) L DDE 
model, [b] the modified LDDE model, [c] and [d] the Wolf et all 120031) PLE 
and PDE model, respectively, of the quasar luminosity function. 



tion rate and, assuming that only the quasar and star-forming 
galaxy populations contribute to FhiCz), it can be expressed as 
the sum of these two different contributions, thus 

/i(3-ab) ^ .fe / 

where e,y.q(z/o,z) and e,y.g(i'o,z) represent the quasar and 
galaxy emissivities at the Lyman limit, respectively. In the fol- 
lowing two sections we separately estimate and discuss these 
different contributors. 

8.2. Quasar contribution 

The integrated luminosity density, or emissivity, of quasars 
is defined as 



4>(L,v,z) LdL 



(15) 



where 0(L, i/jZ) is the quasar luminosity function (LF). In the 
literature, the quasar luminosity function has been estimated 
in a variety of rest frame bands and using different analytical 
descriptions. To convert from the observed passbands to the 
flux at the Lyman limit we approximated the quasar spectra 
energy distribution as a double power law (fi, (xv") with the 
foflowing indices: a = -0.3 at 2500-5000 A jMadau et all 
[1999.) and -0.57 at 228 - 2500 A, as estimated in Sect. H 

The quasar luminosity density has to be constrained 
through the results of quasar surveys. However, even though 
the number of known quasars has probably more than ten- 
folded in the last decade, the relevant range in the redshift- 
luminosity plane is still poorly sampled, mainly due to the fact 
that the luminosity density is dominated by the large number 
of faint QSOs which are hard to identify in large numbers. In- 
stead of relying on a single global approximation, we looked 
into the results of different published surveys and their impli- 
cations for the quasar luminosity density, 
(i) IWo lf et al. (2003) presented two different analytical fits 
to the COMBO- 17 quasar luminosity function: A pure den- 
sity evolution (PDE) and a pure luminosity evolution (PLE) 
model, both approximated by third order polynomials. The 
two derived emissivity laws predict a steep decrease towards 



high redshift with very similar slopes, but showing different 
pea k contributions. 

(ii) iBongiorno et al.l ( 120071) used a combined sample of 
quasars from the VIMOS-VLT Deep Survey (VVDS) and 
SDSS, pro bing a larger luminosity interval compared to 
IWolf e t aTl (2003). They employed a luminosity-dependent 
density evolution model (LDDE) with the limitation of a 
fixed bright end slope of the luminosity function. The result- 
ing em issivity is considerably flatter than that of IWolf et al.l 
(l2003|), 

(iii) i Richards et al.l (120061) estimated the luminosity function 
from the third SDSS data release (DR3). Their DR3 quasars 
are so bright that they were able to probe only the bright end 
slope of the luminosity function and they found significant 
evidence for a slope variation over redshift. For our pur- 
poses their luminosity function cannot be used since their an- 
alytical fits cannot be extrapolated to magnitudes fainter than 
Mi < -26. However we included the slope in formation into an 
additional LDDE model obtained from the IBongiorno et al.l 
(i2007i) model, now allowing for a variable bright-end slope. 
The inferred emissivity decline at high redshift is slightly 
steeper, but still not as pronounced as in the two models by 
IWolfet al. (2003). 

Figure [TT] presents a comparison of the emissiviti es esti - 
mated from the lumino sity functions of ,Wolf et aTl (l2003h : 
IBongiorno etaTI (120071) and IRichards et al.1 (120061) . We as- 
sume that these four analytic descriptions bracket the true evo- 
lution of the emissivity. Using Equation[T3]we converted the 
luminosity density into an QSO contribution to Fhi which we 
denote as Fei.q. 

8.3. Contribution by star-forming galaxies 

The second source population of the UV background to 
be considered are galaxies actively forming stars. The de- 
tection of high redshift galaxies is still extremely challenging 
and only a few deep surveys are available for this calculation. 
Among those reaching redshifts larger than z = 3, we consider 
the results by T resse et al. (2007) using the VVDS sample and 
bv lOuchi etaTI (12004) using the Subaru Deep Field (SDF) sur- 
vey. In order to compute the galactic emissivity at the Lyman 
limit starting from the observed UV luminosity density we 
evaluated 

ei.,g(i^,2) = /esc g(i^) ei]v(i^,z). (16) 

where /esc is the escape fraction of UV photons from galaxies, 
euv(i^,z) is the observed galactic UV emissivity per unit co- 
moving volume, and giv) is the Lyman continuum frequency 
dependence of g alaxies. We assumed that gii^) oc v'^ (e.g. 
Ilnoue et ani2006}) at all redshifts and a constant escape frac- 
tion of 10%. These quantities are highly uncertain and still 
controversially debated. However, we note that these assump- 
tions affect the global scaling of the points and not their rela- 
tive positions, unless one introduces a redshift-dependent es- 
cape fraction (which is possible, but adds another poorly con- 
strained degree of freedom). We then used Eq. [13] to con- 
vert the observed emissivities into a stellar contribution to Fhi 
which we denote as FHi.g. Notice that we refrained from at- 
tempting to construct a smooth interpolation (or even extrap- 
olation) formula for eu.g(y,z)- We only evaluated the con- 
tribution by galaxies rHi,g at the points given by the above 
mentioned surveys. 

8.4. Comparison to the observed photoionization rate 
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Figure 13. Non-QSO contribution to the UV background photoionization 
rate. The different dashed Hues represent the required photoionization rates 
which added to the QSO contribution lead to the observed constant photoion- 
ization rate. The solid triangles and circles sho w Fhi estim ated from the 
stellar emissivities of Tresse et al. (2007) and Ouc hi et alj 12004 ) respectively 
and assuming a constant escape fraction of 10%. 



Figure 14. Non -QSO contribution to the UV background photoionization 
rate as in Fig. 1131 The different dashed lines represent the required photoion- 
ization rates which added to the QSO contribution lead to a mildly decline 
in Fhi towards high redshift. The solid triangles and circles show Fm es- 
timated from the stellar emissivities of Tresse et al. (2007) and Ou chi et alj 
<2004i) respectively and assuming a constant escape fraction of 10%. 



We now consider how our measurements of the overall 
photoionization rate compare with the different contributions. 
Figure [12] shows the different models of the quasar contri- 
bution to r along with the estimates based on the proximity 
effect. The quasar contribution shows a peak around z — 2 
and strongly declines towards higher redshifts, however with 
slopes varying between the different models. 

A number of factors make the QSO contribution to the pho- 
toionization rate uncertain. The analytical description of the 
quasar LF is by far the most prominent one, but also the as- 
sumed quasar SED, the slope of the UVB or the mean free 
path of ionizing photons come into play especially in the 
global scaling of Fhi- Despite the differences in the predicted 
rHi.q(z), all models show a substantial decline in the quasar 
contribution to the global photoionization rate, implying that 
galaxies must provide a substantial additional contribution in 
order to maintain the IGM at the observed high ionization 
level. 

In order to compare the UVB with the predictions from 
source counts, we subtracted the QSO contribution from the 
observed FhiCz) obtained from the proximity effect analysis. 
For our SDSS measurements, the UVB is approximately con- 
stant at Fhi = -1 1 .78. The residuals after subtracting the dif- 
ferent luminosity densities are show in Fig. [13] In all cases 
rHi.g(z) is dominating the global photoionization rate for red- 
shifts larger than three, independently of the details of the as- 
sumed QSO contributions (models [a] to [d]). At smaller red- 
shift the QSO c ontribution inferred from the PLE model by 
I Wolf et all (12003 ) (model [c]) is higher than that of the other 
models, and consequently the stellar contribution is smaller 
for this case. 

We also considered the extreme case of a UVB completely 
dominated by quasars at z 2, by rescaling all QSO emissiv- 
ity models to our proximity effect measurement at z ~ 2. This 
gave us the minimum stellar contribution as a function of red- 
shift. This is zero by design at z = 2, but increases rapidly and 
is in fact very similar to the model [c] in Fig. [13] 

As the uncertainties on the global UVB photoionization rate 
are still considerable, we explored the impact of this uncer- 
tainty by alternatively adopting a slowly declining UVB as 
favored by our UVES results in Paper 11. We computed a lin- 
ear fit to rHi(z) and simply extrapolated this towards higher 



redshift. We then again subtracted the predicted QSO con- 
tribution and plotted the residuals in Fig. [14] This has a sig- 
nificant impact on the predicted stellar contribution. Recall, 
however, that the UVES data from Paper 11 leave the UVB at 
z > 4 essentially unconstrained. 

Fig.[T3]and Fig. [reconstitute two independent assessments 
of the possible stellar contribution to the UVB. If Fhi at 
z ~ 4.2 is as high as suggested by our highest-redshift SDSS 
data point, it would be fully consistent with the predictions 
based on the observed UV luminosity density of galaxies 
at 3 < z < 4. Extrapolating our nearly flat UVB towards 
even higher redshifts would then result in an increasing UVB 
deficit based on the observed galaxy surveys. However, that 
extrapolation is unsupported by the coverage of our data. Fur- 
thermore, some degree of deficit is also perfectly conceivable 
since high-redshift surveys as the SDF sample only the lu- 
minous part (considerably above L*) of the LBG population, 
and the number density of smaller star forming galaxies is still 
very poorly known. 

If on the other hand the UVB should continue to decline 
towards higher z as obtained by extrapolating the results from 
our UVES sample, the counts of star forming galaxies are in 
good agreement at z > 4, but almost too high in the better 
constrained range of 3 < z < 4. We reiterate that these conclu- 
sions rest on the assumption of a constant UV escape fraction. 
By adjusting /esc and allowing for a variation with z, the stellar 
contribution to the UVB becomes essentially unconstrained. 

9. CONCLUSIONS 

Quantifying the metagalactic Hi photoionization rate at 
z > 4 remains a challenge. The three available methods all 
provide estimates that are not inconsistent with each other, 
but the uncertainties are still considerable. For example, the 
poorly constrained escape fraction of UV photons in galaxies 
makes a safe prediction of the UVB from population model- 
ing virtually impossible at present. Likewise, the fundamen- 
tally very powerful method to match numerical simulations 
to the observed opacity evolution of the UVB has to rely on 
several theoretical assumptions that influence the outcome of 
the simulations. The proximity effect is arguably the most 
direct way to measure a quantity that can be translated into 
an estimate of Fhi. The fundamental issue here is the ques- 
tion whether large-scale H l overdensities around quasars can 
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weaken the proximity effect and thus bias the UVB towards 
high values. Our SDSS sample provides clear evidence for the 
presence of such overdensities, but our analysis also suggests 
that the fraction of QSOs heavily affected by overdensities is 
small. Certainly our data point at z ~ 4.2 is very robust and 
unlikely to be affected by any significant overdensity bias. 

Compared to previous analyses of the proximity effect in 
quasar spectra, the main improvement of the present study 
lies the size of the dataset. It constitutes an increase by a 
factor of ~ 20 compared to the hitherto largest sample by 
IScott et aP (2000), at similar spectral resolution, and by a 
factor of 50 compared to our recent evaluation of 40 high- 
resolution spectra from VLT/UVES (Paper II). At redshifts 
z > 4, our new SDSS dataset contains 34 objects, whereas pre- 
vious investigations of the proximity effect at such high red- 
shifts were always based on individual obj ects (IWilliger et alJ 
[Tool: iLu et alJ[!996l: iSavaglio et al!'T997fl. 

We have also introduced an automated continuum deter- 
mination and correction algorithm, which is of fundamental 
importance in quantifying the inevitable effects of continuum 
misplacement. The quality of this correction is impressively 
demonstrated by the fact that our inferred global evolution 
of the Lya forest optical depth is very close to the results of 
most recent studies using high spectral resolution. We thus 
conclude that the moderate spectral resolution of the SDSS 
data has been well accounted for and is not an issue. 

In continuation of our previous work on the proximity ef- 
fect (Papers I and II) we have searched for the proximity effect 
signature not only in the combined sample, but in individual 
quasar spectra. We detected the effect in 98 % of the objects, 
with the remaining 2 % all clearly affected by strong asso- 
ciated absorption. With a sample of this size we could con- 
struct the 'proximity effect strength distribution' (PESD) in 
several redshift bins, with 30-600 QSOs contributing to each 
bin. The observed proximity effect strengths display a consid- 
erable dispersion and a pronounced asymmetry. Both the dis- 
persion and the skewness of the PESD decrease with increas- 
ing redshift, just as predicted from a simple model of ran- 
domly distributed absorption lines following the global line 
density redshift evolution. While it is clear that quasar ab- 
sorption lines are not exactly randomly distributed but follow 
the cosmic web, the distribution of pixel optical depths can be 
reasonably well described as a random process. In Paper II we 
predicted a very strong change of the PESD width with red- 
shift, but owing to a sample size of only 40 objects we could 
not demonstrate the reality of the effect. The SDSS data show 
this trend very clearly. 

The fact that the PESD at higher redshift becomes increas- 
ingly narrow has profound implications for the accuracy of 
measuring the UV background. While any QSO at z ^ 2 may 
show a proximity effect strength that deviates by up to 2 or- 
ders of magnitude from the expectation value, the dispersion 
is much lower at z > 4, and relatively small samples, even in- 
dividual objects can give meaningful constraints on the UVB. 
This is impressively confirmed by the subset of 34 QSOs at 
z > 4 in our sample. 

While our measurement of the metagalactic H I photoion- 
ization rate at z ^ 4 has yielded a rather high value, it is never- 
theless fully consistent with most other recent estimates, cer- 
tainly within 1-2(7 deviations. There is thus a clear trend that 
the UVB is almost constant at log Thi = -1 1 .78 ± 0.07 s"' be- 
tween z — 2 and z — 4, at best it decreases very mildly. It 
will be a most challenging task in the future to push the cur- 



rent limit of z — 4 towards higher redshifts and see whether 
the UVB continues to be flat, or whether a downturn will be 
detected. 
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